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Sector decomposition in its practical aspect is a constructive method used to evaluate Feynman 
integrals numerically. We present a new program performing the sector decomposition and in- 
tegrating the expression afterwards. The program can be also used in order to expand Feynman 
integrals automatically in limits of momenta and masses with the use of sector decompositions 
and Mellin-Barnes representations. The program is parallelizable on modern multicore com- 
puters and even on multiple computers. Also we demonstrate some new numerical results for 
four-loop massless propagator master integrals. 
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1. Introduction 

Originally sector decomposition was used as a tool for analyzing the convergence and prov- 
ing theorems on renormalization and asymptotic expansions of Feynman integrals [[[], |2|, ||, f|, Bj|. 
After [0], the sector decomposition approach has become an efficient tool for numerical evaluating 
Feynman integrals (see Ref. ^ for a recent review). At present, there are two public codes per- 
forming the sector decomposition []8|] and [Q]. The latter one was developed by the present authors; 
it is named FIESTA which stands for "Feynman Integral Evaluation by a Sector decomposiTion 
Approach". Recently FIESTA was greatly improved in various aspects Jlo|]. 

During the last year FIESTA was widely used, some of application are listed in [juj]. In [12] 



we used FIESTA in order to confirm numerically the recent analytic results for master integrals 



(Mi's) for four-loop massless propagators which recently were analytically evaluated in [ 13]. Here 



we provide some more numerical results for extra orders in epsilon expansions for these Mi's. 

2. Theoretical background and software structure 

FIESTA calculates Feynman integrals with the sector decomposition approach. After per- 
forming Dirac and Lorentz algebra one is left with a scalar dimensionally regularized Feynman 



A A 

integral [14] F(a\, ... ,a n ) J ■■■ J ,),'••• „„ ' , where d = 4 — 2e is the space-time dimension, a n are 

fct ...b* n 



indices, / is the number of loops and \/E n are propagators. We work in Minkowski space where 
the standard propagators are the form 1 / (m 2 — p 2 — iO). Other propagators are permitted, see [0]. 
Substituting = daa a '~ l e~' Eia , after usual tricks [^J, performing the decomposition of 

the integration region into the so-called primary sectors [§] and making a variable replacement, 
one results in a linear combination of integrals J^. =0 dxi . . .dx n > ^ITj=i x /' J ^pl'-id/J' 

If the functions na d no singularities in e, one would be able to perform the expansion 

in £ and perform the numerical integration afterwards. However, in general one has to resolve 
the singularities first. Thus, one starts a process the sector decomposition aiming to end with a 
sum of similar expressions, but with new functions U and F which have no singularities (all the 
singularities are now due to the part ]Tj=i* / /' )• The way sector decomposition is performed is 
called a sector decomposition strategy ([||, [| |9]]) and is an essential part of the algorithm (let us 



also mention a geometrical approach to sector decomposition [15] which is rather complicated in 
implementation as a strategy on a computer but promises to be the optimal one). 

After the sector decomposition one resolves the singularities by evaluating the first terms of 
the Taylor series: in those terms one integration is taken analytically. Afterwards the £-expansion 
can be performed and finally one can do the numerical integration. 

FIESTA is written in Mathematica [ |l7| ] and C. The user is not supposed to use the C 
part directly as it is launched from Mathematica via the Mathlink protocol in order to per- 
form a numerical integration. To run FIESTA, the user has to load the FIESTA package into 
Mathematica 6 or 7. In order to evaluate a Feynman integral one has to use the command 
SDEvaluate [UF [ loop_momenta, propagators, subst] , indices, order] , where 
loop_momenta is a list of all loop momenta, propagators is a list of all propagators, subst 
is a list of substitutions for external momenta, masses and other values. For example, 
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SDEvaluate [UF [ {k} , { -k 2 , - (k+pi) 2 , - (k+pi+p 2 ) 1 , - (k+pi+p 2 +P4) 2 } , 
{p 2 ->0, p 2 ->0, P4 ->0, pi p 2 ^-s/2,p 2 P4-)— t/2,pi P4 (s+t) /2, 
s-)--3,t-»-l}], {1,1,1,1},0] 
evaluats the massless on-shell box diagram with Mandelstam variables equal to —3 and — 1. 



3. Numerical results for four-loop massless propagators 
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Figure 1: M^i-M^: the thirteen complicated four-loop master integrals according to [pj[]. The two Mi's 
M52 and M43 can be identically expressed through the three-loop nonplanar MI A^>- 



In [18] a full set of four-loop massless propagator-like MTs was identified. There are 28 inde- 



pendent MTs. Analytical results for these integrals were obtained in J13Q. The most complicated 
Mi's are demonstrated on Fig. |l|. e m after My stands for the maximal term in £-expansion of My 
which one needs to know for evaluation of the contribution of the integral to the final result for a 
four-loop integral after reduction is done, see [|TJ|]. Two of the complicated integrals (M43 and M52) 
are related by a simple factor with the three-loop MI [ |l2"l ] so it is enough to evaluate remaining 
eleven complicated Mi's M^-M^ as well as first three terms of the £-expansion of Nq. 

We calculated them (for q 2 = —1) using FIESTA with the Cuba[16] Vegas integrator and 
1 500 000 sampling points for integration. Our results alongside with the corresponding analytical 
expressions (transformed to the numerical form) from [|13|] look like follows 1 : 

M 34 £~ 4 : 0.08333 ± (0.08333); £~ 3 : 0.916667 ± 0.000018 (0.91666); £~ 2 : 5.64251 ± 0.00022 
(5.6425109); e -1 : 27.6413 ± 0.00077 (27.6412581); £°: 98.638 ± 0.0034 (98.637928); e 1 : 
342.736 ± 0.012 (342.7349920); £ 2 : 857.88 ± 0.048 (857.8735165); £ 3 : 2659.84 ± 0.19 
(2659.825402); £ 4 : 4344.28 ± 0.75 (unknown); £ 5 : 17483.1 ± 5.7 (unknown). 



M 35 £- 2 : 0.601028 ± 0.000012 (0.601028); e _1 : 7.4231 ± 0.00024 (7.423055); £ u : 44.9127 ± 
0.00073 (44.91255); e 1 : 217.023 ± 0.0037 (217.0209); £ 2 : 780.436 ± 0.013 (780.432); £ 3 : 
2678.13 ± 0.053 (unknown); £ 4 : 7195.9 ± 0.3 (unknown). 



'Please, note that the overall normalization used by FIESTA is different from the one employed by the authors of 
@, see 
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M 36 e _1 : 5.184645 ± 0.000042 (5.1846388); e°: 38.8948 ± 0.00039 (38.8946741); e 1 : 240.069 
± 0.0019 (240.0684359); e 2 : 948.623 ± 0.0091 (unknown); £ 3 : 3679.77± 0.06 (unknown). 

Mm e- { : 20.73860 ± 0.00023(20.7385551); £°: 102.033 ± 0.003 (102.0326759); e 1 : 761.60 ± 
0.011 (761.5969858); £ 2 : 2326.18 ± 0.062 (unknown); £ 3 : 12273.6 ± 0.4 (unknown). 

M 42 e _1 : 20.73860 ± 0.00024 (20.7385551); £°: 145.381 ± 0.0029 (145.3808999); e 1 : 985.91 
± 0.014 (985.9082306); £ 2 : 3930.65 ± 0.076 (unknown); £ 3 : 17486.6 ± 0.6 (unknown). 

M 44 £°: 55.58537 ± 0.00031 (55.5852539); e 1 : 175.325 ± 0.004 (unknown); £ 2 : 1496.52 ± 0.02 

M 45 £°: 52.0181 ± 0.0003 (52.0178687); e 1 : 175.50 ± 0.0036 (175.496447); £ 2 : 1475.272 ± 
0.0098 (unknown); £ 3 : 2623.5 ± 0.1 (unknown). 

M S i e _1 : -5.184651 ± 0.000048 (-5.184638); e°: -32.0962 ± 0.00057 (-32.09614); e 1 : -91.158 
± 0.0052 (-91.1614); £ 2 : 119.06 ± 0.043 (unknown); £ 3 : 2768.6 ± 0.45 (unknown). 

Af £°: 20.73857 ± 0.00026 (20.7385551); e l : 190.60 ± 0.0023 (190.600238); £ 2 : 1049.20 ± 
0.014 (1049.194196); £ 3 : 4423.84 ± 0.072 (unknown); £ 4 : 16028.8 ± 0.5 (unknown). 

M 61 e _1 : -10.36931 ± 0.00006 (-10.3692776); £°: -70.990 ± 0.0011 (-70.99081719); e 1 : - 
21.650 ± 0.013 (-21.663005); £ 2 : 2832.69 ± 0.096 (unknown). 

M 62 e -1 : -10.36933 ± 0.00006 (-10.36927); £°: -58.6187 ± 0.0013(-58.6210); e 1 : 244.681 ± 
0.015 (unknown). 

M 63 e- 1 : -5.18467 ± 0.000042 (-5.184638); £°: 14.3989 ± 0.00081 (14.39739); e 1 : 739.979 ± 
0.0099 (unknown). 



Here for each MI we provide our numerical result for coefficients of £-expansion in comparison (in 
parentheses) with the known from [|l3|] analiycal results (if any). As we can see, our calculations 
reproduce the result of [13] with 3-4 correct digits. The extra terms in the £-expansion of each MI 
which are currently unavailable analytically but are necessary for future five-loop calculations. 



4. Conclusion 

Usually, analytical evaluation of multiloop MI is a kind of art. It requires a lot of efforts (and 
CPU time). In many situations, independent checkup is hardly any possible in reasonable time. 
That is why the simple in use tools for numerical evaluation like FIESTA are important. 

Acknowledgments. This work was supported in part by DFG through SBF/TR 9 and the Rus- 
sian Foundation for Basic Research through grant 08-02-01451. 
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